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ABSTRACT 

We report the confirniation and characterization of a transiting gas giant planet orbiting the M 
dwarf KOI-254 every 2.455239 days, which was originahy discovered by the Kepler mission. We use 
radial velocity measurements, adaptive optics imaging and near infrared spectroscopy to confirm the 
planetary nature of the transit events. KOT254b is the first hot Jupiter discovered around an M-type 
dwarf star. We also present a new model-independent method of using broadband photometry to 
estimate the mass and metallicity of an M dwarf without relying on a direct distance measurement. 
Included in this methodology is a new photometric metallicity calibration based on J — K colors. We 
use this technique to measure the physical properties of KOT254 and its planet. We measure a planet 
mass of Mp sin i = 0.505 Mjup, radius Rp = 0.96 Rjup and semimajor axis a = 0.030 AU, based on 
our measured stellar mass Mi, — 0.59 and radius Ri, = 0.55 R©. We also find that the host star 
is metal-rich, which is consistent with the sample of M-type stars known to harbor giant planets. 



1. INTRODUCTION 

Of the hundreds of planets discovered outside of the 
solar system the class known as hot Jupiters, with peri- 
ods P < 10 days and Mp > 0.2 Mjup, has provided much 
of our knowledge about the internal structures and com- 
positions of exoplanets (Charbonneau et al. 2007; Winn 
2008). The crucial role hot Jupiters play in present-day 
exoplanetary science is not because they are common; 
the occurrence rate of hot Jupiters in the solar neigh- 
borhood is only about 1% (Marcy et al. 2005). Rather, 
close-in giant planets are the most easily detectable vari- 
ety for Doppler and transit surveys owing to the large re- 
flex motions and transit depths they induce. Once found 
to transit their host stars, their relatively large atmo- 
spheric cross-sections make them favorable targets for 
various follow-up studies to measure planet properties 
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such as atmospheric composition, atmospheric temper- 
ature profiles, phase curves, and spin-orbit angles (e.g. 
Knutson et al. 2009; Crossfield et al. 2010; Madhusud- 
han & Seager 2011; Albrecht et al. 2011). 

While only one in every hundred Sun-like stars harbors 
a close-in gas giant planet, the occurrence of hot Jupiters 
is even further depressed around the Galaxy's most nu- 
merous denizens, the M dwarfs. This empirical finding 
has emerged from various Doppler surveys of M dwarfs, 
which have detected zero close-in giant planets among 
roughly 300 target stars with masses Mi, < 0.6 Mq, de- 
spite the ready detectability of the large Doppler am- 
plitudes of these short-period giants (Endl et al. 2003; 
Johnson et al. 2010a). This lack of hot Jupiters around 
M dwarfs, as well as the overall dearth of giant planets 
around low-mass stars, is likely due to the inefficiency 
of the planet-formation process within low-mass proto- 
planetary disks (Laughlin et al. 2004; Ida & Lin 2005; 
Kennedy & Kenyon 2008). This notion is further bol- 
stered by the elevated occurrence of giant planets around 
stars more massive than the Sun (Johnson et al. 2007, 
2010a). 

To date, no planet more massive than 0.1 Mjup has 
been discovered around an M dwarf with a period less 
than 30 days^'^. However, orbital companions at either 
mass extreme have been detected around M dwarfs, and 
some have been found to transit their stars. At the high- 
mass end, two transiting brown dwarfs have been found 
by the MEarth and Kepler transit surveys, respectively 
(Irwin et al. 2010; Johnson et al. 2011a). At the other 
end of the mass scale, transits of the RV-detected, short- 
period Neptune Gl 436 b were detected by ground-based 
follow-up photometry (Butler et al. 2004; Gillon et al. 
2007), and the MEarth survey discovered a transiting 
super earth orbiting the nearby, low-mass star GJ 1214 
(Charbonneau et al. 2009). However, no close-in planet 
with a mass intermediate to those previous discoveries 
has been found around an M dwarf. This leaves a gap 

http://exoplanets.org 
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Fig. 1. — Normalized Kepler photometry vs. time showing the 
rotational modulation of star spots, from which we derive a rotation 
period. 

in the companion-mass continuum that spans more than 
two orders of magnitude, from 0.1 Mjup to 30 Mjup. 

One of the planet candidates discovered by the Kepler 
mission, KOI-254. 01, provides an opportunity to fill in 
this missing region of parameter space. Borucki et al. 
(2011) reported a transit depth oi 6 — 3.9093% and a 
period of P = 2.455239 days for the planet candidate. 
The large transit depth and lack of a visible secondary 
eclipse, together with the lack of a pixel shift in the stellar 
position during transit, reduces the a priori probability 
that the transit is a false-positive to 1.9% (Morton & 
Johnson 2011). However, it is nonetheless desirable and 
prudent to acquire a full suite of follow-up observations 
to rule out a false positive. We also characterize the 
stellar characteristics of the low-mass host star in order 
to accurately and precisely measure the properties of the 
planet. 

This contribution thus builds on the work of Muirhead 
et al. (2011), hereafter Paper 1, who presented a spec- 
troscopic analysis of all of the M dwarf KOIs announced 
by Borucki et al. (2011). We add to that work by in- 
cluding additional data from our multi-band, multi-site 
follow-up campaign for the specific case of the KOI-254 
planetary system. In § 2 we present our observations of 
KOI-254. In § 3 we describe our analysis of the Kepler 
light curves, KIC photometry, infrared spectra, and ra- 
dial velocity orbit. We describe our results in § 4 and 
conclude with a brief discussion and summary in § 5. 

2. OBSERVATIONS 
2.1. Kepler Photometry 

The Kepler Object of Interest KOI-254 
(=KIC 5794240, =2MASS 19312949-^4103513) is 
identified in the Kepler Input Catalog (KIC; Batalha 
et al. 2010; Brown et al. 2011) as a cool dwarf star with 
Tcff — 3948 K and log 5 = 4.54 and a Kepler magnitude 
Kp = 15.98 (no errors are reported). Periodic, deep 
transit events were detected and reported by Borucki 
et al. (2011), with a period P = 2.455239±0.000004 days 
and a depth of 3.9093%. 

KOI-254 was not observed during the QO observing 
quarter, but was during subsequent quarters. We down- 
loaded the Ql and Q2 data sets from the NASA Mul- 
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timission Archive at STScI (MAST). The photometric 
measurements span 2009 May 13 to 2009 Sept 16 and 
have a 29.4244-minute cadence and a median fractional 
uncertainty of 6 x 10~^. In addition to the transit events, 
the data show a clear 14-day periodicity with peak-to- 
peak variations of approximately 5% persisting coher- 
ently over the entire 122-day time baseline. This addi- 
tional variability is most likely due to rotational modu- 
lation of star spots, which is useful for setting limits on 
the system age (§ 4.3). 

The KIC lists additional photometric data in the Sloan 
Digital Sky Survey (SDSS) bands, as well as near infrared 
(NIR) photometry from the Two-Micron All-Sky Survey 
(2MASS Skrutskie ct al. 2006), which we reproduce in 
Table 1. 

2.2. Nickel Z-band Photometry 

To check for achromatic transit depths and to update 
the ephemeris over that provided by Q2 data, we ob- 
served the transit predicted to occur on UT 2011 June 
30 using the 1-m Nickel telescope at Lick Observatory 
on Mt. Hamilton, California. We used the Nickel Direct 
Imaging Camera, which comprises a thinned Loral 2048^- 
pixel CCD with a 6'.3 square field of view (e.g., Johnson 
et al. 2008). We observed through a Gunn Z filter, used 
2x2 binning for an effective pixel scale of 0'.'37 pixel" ^, 
and a constant exposure time of 120 seconds. We used 
the fast readout mode, with approximately 13 s between 
exposures to read the full frame and reset the detector. 

The conditions were clear, and the seeing and telescope 
optics delivered a full width at half maximum of ^ 1'.'5. 
We began observing as soon as possible after sunset at an 
airmass of 1.55 and observed continuously for 7.4 hours 
bracketing the predicted transit midpoint, ending at an 
airmass of 1.21. Guiding was unsteady for the first 1.5 
hours but settled down thereafter, subsequently needing 
only small occasional adjustments to keep the stars as 
nearly as possible on the same pixels. 

We measured the instrumental magnitude of KOI-254 
with respect to four nearby stars with KIC numbers 
5794268, 5794279, 5794302, and 5794355. An aperture 
width of 22 pixels (with a sky annulus of inner and outer 
diameter of 27 and 32 pixels, respectively) gave the lowest 
out-of-transit photometric scatter. The Nickel dome par- 
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Fig. 2. — Keck/HIRES radial velocity measurements of KOI-254 
(black with error bars) phased to the transit-based orbital period. 
The error bars represent the quadrature sum of the internal mea- 
surement errors and 23 m s~^ of jitter based on our analysis of 
low signal-to-noise measurements of a stable M dwarf. The dashed 
line shows the best-fitting Keplerian model, and the red squares 
are the model evaluated at the times of measurement. 

tially occults the telescope when observing within 5.7 de- 
grees of zenith; this introduced systematic photometric 
variations that caused us to excise 40 minutes of obser- 
vations, including the transit egress. We converted the 
Nickel timestamps to BJDutc using the techniques of 
Eastman et al. (2010) to be consistent with the Kepler 
data. The Nickel photometric measurements, phased at 
the Kepler transit period, are shown in Figure 7. 

2.3. Keck/HIRES Velocimetry 

We obtained spectroscopic observations of KOI-254 
at Keck Observatory using the High-Resolution Echelle 
Spectrometer (HIRES) with the standard iodine-cell 
setup used by the California Planet Survey (Howard et al. 
2010). Because of the star's faintness {V = 16.88) we 
used the C2 decker corresponding to a projected size 
of 14'.'0x0'.'851 to allow sky subtraction and a resolving 
power of i? = A/AA « 55,000. We obtained 14 obser- 
vations of KOI-254, all with an exposure time of 1200 
seconds, resulting in a signal-to-noise ratio (SNR) of 13- 
16 at 5500 A. 

We used the Butler et al. (1996) iodine cell method 
to measure the radial velocity time series of KOI-254 
with respect to an iodine-free "template" observation 
that has had its instrumental profile removed through 
deconvolution. However, since KOI-254 is so faint, ob- 
taining a high-SNR, high-resolution template would be 
prohibitively expensive in observing time. We instead 
selected a surrogate star with the same spectral type as 
KOI-254 (HD 199305, MOV) and adjusted the line depths 
and widths using the "spectral morphing" technique of 
Johnson et al. (2006). We then used this morphed, de- 
convolved template in our RV analysis. 

To ensure that our Doppler analysis pipeline produced 
reliable RV measurements, we obtained a sequence of 
test exposures of the MOV star HIP 36834, which has 
demonstrated RVs stable at the 5.3 m s^^ level over the 
past seven years at Keck. Our test exposures were made 
at SNR =15 and a morphed template based on the same 
surrogate stellar spectrum used for KOI-254. We find a 
root- mean square (rms) of 23 m s~^ during a time span 
of 20 days. We adopt this rms as our RV measurement 
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^ Uncertainties represent formal mea- 
surement uncertainties and do not in- 
clude the 23 m s~^ instrumental sys- 
tematics we include in our orbit analy- 
sis. 

uncertainty for KOI-254. 

Figure 2 shows our 14 RV measurements of KOI-254 
phased at the transit period, with the error bars corre- 
sponding to our estimated 23m measurement pre- 
cision. The best-fitting Keplerian model is shown as a 
dashed line and described in § 3. The RV measurements, 
HJD times of observation and internal uncertainties are 
listed in Table 2. 

2.4. Keck Adaptive Optics Imaging 

We obtained adaptive optics (AO) images of KOI-254 
on June 24, 2011 UT using the NIRC2 camera at Keck 
II, in order to reduce the likelihood for false positives 
by searching for any sources that could mimic a plan- 
etary transit signal, such as a nearby eclipsing binary. 
With r = 16.11 mag, KOI-254 is relatively faint for nat- 
ural guide star observations. Nevertheless, the conditions 
were excellent with very little cirrus and seeing « 0'.'5. 
We were able to close the AO system control loops on the 
star with a frame rate of 30 Hz. With sufficient counts 
in each wavefront sensor subaperture, a stable lock was 
maintained for the duration of the observations. 

We acquired a sequence of 9 dithered images in the K' 
filter (central A = 2.12/zm) using the NIRC2 medium 
camera (plate scale=20 maspix"-'^). Each frame con- 
sisted of 6 coadds with 10 s of integration time per 
coadd, totaling 9 minutes of on-source exposure time. 
The images were processed by removing hot pixels, sub- 
tracting the sky-background, and aligning and coadding 
the cleaned frames. Figure 3 shows the final reduced 
AO image. The field of view is 9'.'8 x 9'.'8, which cor- 
responds to 2.45 Kepler pixels on a side. No obvious 
contaminants were identified in the immediate vicinity 
of KOI-254. Using the limits from our AO images, we 
are able to rule out any contaminating sources to a level 
of AK' = {1.0,3.2,5.6,6.9,6.9} at 3ct for separations of 
{0'.'25, 0'.'5, I'.'O, 2"0, 4'.'0} respectively 

2.5. Palomar Near Infrared Spectroscopy 

We obtained near-infrared spectra of KOI-254 with 
the TripleSpec Spectrograph at the Palomar Observatory 
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Fig. 3. — Contrast limits based on the high spatial resolution K' 
image of KOI 254 taken with NIRC2 and the Keck AO system on 
June 24th, 2011 UT. The contrast is measured in magnitudes with 
respect to the central star. The AO observations probe regions 
interior to the Kepler photometric aperture and help to eliminate 
candidate false-positive signals at brightness levels comparable to 
the transit depth. No nearby stars were detected. 

200-inch Hale Telescope (Herter et al. 2008). The spec- 
tra were obtained as part of a survey of low-mass Kepler 
Objects of Interests described in Paper 1. TripleSpec is 
a cross-dispersed, long-slit, ncar-infrarcd spectrograph, 
dispersing a 1 x 30 arcsecond slit from 1.0 to 2.5 /im 
across 5 orders at resolution A/ A A = 2700. Two posi- 
tions on the slit, A and B, were used for each target, 
and exposures were taken in an ABBA pattern. Obser- 
vations were taken on June 12, 13 and 18 of 2011. On 
June 12th, the seeing was photometric, and 3 ABBA sets 
were taken with 30 exposures at each position. On June 
13th, seeing was 1.5 to 2 arcscconds, and 4 ABBA sets 
were taken with 60 exposures at each position. On June 
18th, seeing was back to photometric, and 4 ABBA sets 
were taken with 60 sec exposures. 

Spectroscopic observations at near-infrared wave- 
lengths must contend with telluric absorption lines in- 
troduced by the Earth's atmosphere. Telluric lines vary 
strongly with airmass and humidity, and are calibrated 
by observing an object with a known spectrum at sim- 
ilar airmass and near in time to the science target ob- 
servations. To calibrate the telluric lines in KOI-254 we 
observed HD 183204, an AO star found using SIMBAD, 
which was within 0.1 airmasses of KOI-254 at the times 
of the observations. 

The spectra were extracted using a version of the Spex- 
Tool program modified for the Palomar TripleSpec Spec- 
trograph (Gushing et al. 2004). The xtellcor pack- 
age within SpcXTool accepts spectra of AO stars and 
compares them to a model spectrum of Vega to identify 
and remove telluric absorption lines in a target spectrum 
(Vacca et al. 2003). Figure 4 plots the NIR spectrum of 
KOI 254 along with templates of similar spectra type 
(HD 199305 [MOV] and G1229A [MIV]), with the rele- 
vant spectral features indicated. We only show J H and 
K bands, as the regions between are heavily obscured by 
telluric water vapor. The templates are taken from the 
SpeX library (Gushing et al. 2005; Rayner et al. 2009). 
The NIR spectrum of KOI-254 is consistent with an MO 
dwarf. 

3. ANALYSIS 



Fig. 4. — TripleSpec NIR spectra in Ks band for KOI-254 and 
two other early M dwarfs for comparison (HD 199305 [MOV] and 
G1229A [MIV]). The comparison spectra are taken from the Infra- 
Red Telescope Facility (IRTF) spectral library {gray, Gushing 
et al. 2005; Rayner et al. 2009). We also indicate strong atomic 
features identified by Gushing et al. (2005). if-band contains sig- 
nificantly fewer atomic features than J-band, though the ii'-band 
Na I doublet and Ga I triplet have been shown to be sensitive to 
metallicity (Rojas-Ayala et al. 2010). 

3.1. Analysis of Broadband Photometry 

The process of measuring the physical properties of 
stars generally involves comparing an observed tempera- 
ture indicator (e.g. color) and luminosity (apparent mag- 
nitude, bolometric correction and parallax) to theoretical 
models of stellar evolution (e.g. Nordstrom et al. 2004; 
Valenti & Fischer 2005; Takeda et al. 2008). Such a com- 
parison is typically performed through an interpolation 
of the observables onto a grid of models, and locating the 
combination of stellar properties that provide the clos- 
est match to the observations. However, the process is 
complicated for low-mass stars in general because of in- 
accuracies in theoretical evolution models (Torres 2007). 

Fortunately, for ages > 100 Myr and spectral types ear- 
lier than M5, M dwarfs reside very close to the zero- age 
main sequence, and at a fixed metallicity have a one-to- 
one mapping between mass and luminosity. Indeed, in 
the NIR J, H and Ks bands, the mass-luminosity rela- 
tionship is independent of metallicity (Allard et al. 1997; 
Delfosse et al. 2000). The metallicity dependence of the 
mass-luminosity relationship in optical bands provides 
a means of estimating stellar metallicitics using broad- 
band photometry (Bonfils et al. 2005; Johnson et al. 
2009; Schlaufman & Laughlin 2010). Additional relation- 
ships have been discovered between optical-NIR colors, 
and spectral type (West et al. 2005). These various re- 
lationships each provide leverage in the determination of 
the stellar properties, all without appealing to theoretical 
stellar structure and atmosphere models. 

Our methodology is similar to the approach used by 
Johnson et al. (2011a) to measure the properties of low- 
mass Kepler target star, LHS 6343. In that case Johnson 
et al. took advantage of the common distance and com- 
position of each component of the visual M+M binary 
system to solve for the mass and radius of the binary 
component transited by a brown dwarf. In the present 
work we generalize the methodology for single M dwarfs 
by enlarging the number of color and magnitude relation- 
ships included in the analysis, including a new NIR color- 
based metallicity calibration described in Appendix A. 

We treat the stellar properties {a} = {M^, [Fe/H],(i} 
as free parameters in a set of models that reproduce the 
observables {D} = {fo, -Di, £>jv-i}, where is the 
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star's mass in solar units, [Fc/H] is the metallicity and 
d is tlic distance from the Sun measured in pc. The ob- 
served quantities {D} = {J,H, Ks, V,r— J,aR,r — i} are 
the apparent magnitudes, colors and transit parameters, 
with their associated measurement uncertainties {a}. 

Each observed quantity and its uncertainty are repre- 
sentative of a probability density function, approximated 
by a normal function with mean Di and width ai, which 
can be evaluated at the value predicted by a model as 
a function of the parameters {a}. This calculation pro- 
vides a probability of a datum conditioned on the model 
parameters, and the "best-fitting" values of the param- 
eters are those that maximize the probability of the full 
set of observed data. For the case in which the pdf of 
each measurement is described adequately by a normal 
distribution, one recovers the special situation in which 
maximizing the probability of the data is equivalent to 
minimizing x^, where 
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In this expression, fi{a) is a model function that trans- 
forms the parameters {a} into a prediction of an ob- 
servation Di. The p{aj) are probabiity terms that en- 
code prior knowledge about the parameters {a}. For 
most of these we adopt normal distributions such that 



\lnp{aj)\ = (ap 



dj) /a^., where aprior,j is the most 



likely value of Uj based on prior knowledge, and is 
the width of the normal distribution centered on aprior.j- 
These terms thus serve as "penalty functions" that in- 
crease for deviations far from prior knowledge of the 
parameter values. 

3.2. Relationship Between Aapparent Magnitude and 

{AU,d} 

Delfosse et al. (2000, ; hereafter DOG) provide mass- 
luminosity relationships based on low-mass eclipsing bi- 
naries in the three 2MASS bands {J, H, Ks), with the 
tightest correlation between stellar mass and absolute 
Ks magnitude, Mks- Rather than computing a mass 
based on an absolute magnitude, which we lack due to 
an unknown parallax, we instead evaluate predictions of 
the apparent magnitude for j = {J,H,Ks} given a 
stellar mass M^,, and distance d. This can be expressed 
as 



m,(M^, d) = Mj{M^) + 5 logio(rf/10) + Aj (2) 

where M.j{Mi,) are provided by Jll as polynomial ex- 
pressions with coefficients listed in their Table 3, and Aj 
are extinction terms. We adopt a normally-distributed 
prior distribution for these terms based on Ay = 0.195 ± 
0.02 from the KIC. Further, we adopt the reddening 
law given for sight line to KOI-254 as listed in the 
NASA/IPAC Extragalactic Database based on Schlegel 
et al. (1998). However, we find that our results are un- 
changed if we adopt Ay = 0, indicating that extinction 
is minimal at the distance and direction of KOI-254. 

3.3. Relationship Between Aapparent Magnitude and 
{M,,[Fe/H],4 



The t^-band magnitude can be related to the stellar 
mass and metallicity in a manner similar to the method 
used by Jll, namely 

V{M^, [Fe/H],d)=Mv{M^, [Fe/H]) 

+ 51ogio(d/10) + Ay (3) 

where 



E &v,k[AlKs(M^) + (0.71 [Fe/H] -0.091)]^= (4) 

fc=i 

where the coefficients {&v} are given in Table 3 of Jll. 
The second term in brackets, the linear function of F, 
differs from that used by Jll owing to a revised broad- 
band metallicity calibration that compares well with that 
of Schlaufman & Laughlin (2010). 

An additional constraint on the metallicity is provided 
by om- new J — Ks metallicity calibration described in 
Appendix A and expressed in Eqn. A3. 

3.4. Relationship Between a/ Ri, and M^, 

The mass can also be estimated using the scaled semi- 
major axis from the fit to the transit light curve (e.g. 
Seager & MaUen-Ornelas 2003; Sozzctti et al. 2007; Winn 
2008). The scaled semimajor axis, an = a/R^ is related 
to the stellar density via Kepler's third law: 
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where we have assumed Mp ^ Mi,. The function 
Ri,{Mi,) relates the stellar mass to the radius, with 
both quantities in solar units. Jll used the empiri- 
cal mass-radius relationship of Ribas (2006). However, 
the radii used in this analysis are based on short-period 
eclipsing binaries, for which there is the possibility that 
rotation-induced magnetic and coronal activity, as well 
as metallicity-inhibited convection may lead to anoma- 
lously inflated radii compared to single M dwarfs (Lopez- 
Morales 2007; Kraus et al. 2011). In the present analy- 
sis, we use an empirical mass-radius relationship based 
on the interferometric radii of nearby M dwarfs measured 
by Boyajian et al. (2012, in prep.), which gives a polyno- 
mial relationship that is very similar to the one used by 
Jll based on the eclipsing binaries in the Ribas (2006) 
sample. 

3.5. Relationship Between r — J color and M^, 

Finally, to further constrain the mass of the star, we 
used the color-luminosity relationships of West et al. 
(2005). Specifically, we used columns 2, 5 and 6 from 
their Table 1 to relate Aij to r — J. We then recast the 
relationship in terms of stellar mass using j(M*) from 
Jll to give 



{r-J)[M^]=0.6587xMj{M^)-1.738 + {Ar-Aj) (6) 



5 6 7 8 9 

Absolute K-band Magnitude 

Fig. 5. — The fractional difference between known parallax-based 
distances of nearby M dwarfs, dparallaxi and estimated distances 
from our broadband photometric method, dest- The differences 
are plotted versus the absolute i^g-band magnitude Mk^- The 
rms indicates that our method has a precision of 13.5%, which 
is consistent with the 68.2% confidence region from our MCMC 
analysis of KOI-254. 

3.6. Joint Analyasis of the Stellar Properties 

We searched for the best-fitting vahies of the model 
paranieters {a} and their uncertainties using a Markov 
Chain Monte Carlo (MCMC) algorithm (See, e.g. 
Tegmark et al. 2004; Ford 2005; Winn et al. 2007; John- 
son et al. 2011b). The MCMC technique uses the data 
and priors to explore the shape of the posterior probabil- 
ity density function (pdf ) for each parameter of a model, 
conditioned on the available data. MCMC, particularly 
with the Metropolis-Hastings algorithm, provides an ef- 
ficient means of exploring high-dimensional parameter 
space and mapping out the posterior pdf for each model 
parameter. 

At each "chain link" in our MCMC analysis one pa- 
rameter is selected at random and is altered by drawing 
a random variate from a transition probability distribu- 
tion. In our case we use a normally distributed, pseudo- 
random number as the random variate. If the resulting 
value of the for the trial parameters is greater than 
the previous value, then the set of trial parameters are 
accepted and added to the chain. If not, then the prob- 
ability of adopting the new parameter set is exp(Ax^), 
where Ax^ is the difference in the fitting statistic from 
the previous and current steps. If the current trial is 
rejected then the parameters from the previous step are 
adopted. The width of the transition function determines 
the efficiency of convergence. If it is too narrow then the 
full exploration of parameter space is slow and the chain 
is susceptible local minima; if it is too broad then the 
chain exhibits large jumps and the acceptance rates are 
low. 

The various empirical calibrations we use are not exact 
and have uncertainties owing to the imperfect calibration 
data used to construct them. For example, we refitted 
the Delfosse et al. mass-luminosity relationship using 
their calibration stars and found an rms scatter in the 
absolute magnitudes about the best-fitting polynomial 
relationships. For example, we find that the rms scat- 
ter in Mks from the Delfosse et al. calibration stars is 



0.55 - 
0.53 _ 

200 300 400 -0.4 -0.2 0.0 0.2 0.4 

d [pc] [Fe/H] 

Fig. 6. — Each panel shows the joint, posterior pdfs for two stel- 
lar parameters at a time, marginalized over the remaining param- 
eters. The contours show the iso- probability levels corresponding 
to {68.2, 95, 99.7}% confidence. The best-fitting values are shown 

as solid circles. 

0.18 mag. We adopt this as our uncertainty in Ks for 
KOI-254 in place of the value reported in second column 
of Table 1. Our adopted NIR uncertainties are given in 
column 3 of Table 1 . 

We tested our methodology on a sample of nearby M 
dwarfs with well-measured parallaxes from van Leeuwen 
(2007), 2MASS photometry and SDSS r magnitudes 
from the Carlesberg Meridian Catalog (CMC14 Evans 
et al. 2002). For our test we treated each of these star as 
though their distances, and Ai Ks values were unknown, 
and used the available photometry in our MCMC scheme 
to estimate their distances. The results are shown in 
Figure 5. The fractional, root-mean scatter (rms) of our 
mass estimates about the DOO values is 13.5%, which we 
adopt as our measurement uncertainty. 

3.7. Stellar Properties Measured From Photometry 

We find a stellar mass M^, = 0.59 ± 0.06 Mq and a 
stellar radius i?* = 0.55 ± 0.11 R©, where the value 
in parentheses is the systematic error in our distance 
estimates based on the root-mean-scatter (rms) scatter 
shown in Figure 5. We also find [Fe/H] = -^0.13 ± 0.13 
and d = 333 ± 50 pc. The 68.2% confidence region for 
our MCMC estimate of the distance is consistent with the 
13.5% rms scatter seen in Figure 5. This is an indication 
that our inflated measurement errors on our photomet- 
ric measurements are estimated properly and account for 
the systematic errors in the empirical calibrations we em- 
ploy in our analysis. Figure 6 shows the two-dimensional 
posterior probability distributions for our stellar param- 
eters. 

3.8. Analysis of NIR Spectra 

The moderate-resolution Jf-band spectra of late-type 

stars contain atomic lines highly sensitive to metallicity 
and continuum regions sensitive to effective temperature 
(Covey et al. 2010; Rojas-Ayala et al. 2010). We use 
the spectral indices and calibrations of Rojas-Ayala et 
al. (2011, submitted) to measure the iron abundance 
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[Fe/H], overall metallicity [M/H], and effective tempera- 
ture Teff of KOI-254. The metallicity and effective tem- 
perature relations use the equivalent widths of the Ca I 
triplet and Na I doublet in ii'-band and the H2O-K2 in- 
dex: a measurement of the deformation of the K-hand 
pseudo-continuum using regions dominated by water ab- 
sorption (based on the H2O-K of Covey et al. 2010). 
The metallicity relations are empirically calibrated us- 
ing M dwarfs with wide F-, G- or K-type companions, 
which have metallicity measurements in the SPOCS cata- 
log (Valenti & Fischer 2005). The effective temperatures 
are measured by interpolating the H2O-K2 value and the 
measured overall metallicity [M/H] of a given star onto a 
grid of H2O-K2, [M/H] and T^e calculated using the BT- 
Settl-2010 model spectral spectra of AUard et al. (2010). 
The grid consists of BT-Sctl-2010 spectra for effective 
temperatures between 2200 and 4000 K with 100 K in- 
crements, and [M/H] values of -1.0, -0.5, 0.0, +0.3 and 
-1-0.5. The H2O-K2 is calculated on the model spectra, 
providing well-sampled grid of H2O-K2, [M/H] and Toff 
values. 

We use Monte Carlo simulations to calculate the mea- 
surement uncertainties of the [Fe/H], [M/H] and Teff. We 
calculate the values for 1000 realizations of the spectra, 
each with noise added based on the per-channel error es- 
timates reported by SpeXTool. We take the standard 
deviation in the resultant distribution of [Fe/H], [M/H] 
and Toff as the measurement uncertainty in those values. 
Rojas-Ayala et al. (2011, submitted) estimate a system- 
atic uncertainty of 0.14 for [Fe/H] and 0.10 for '[M/H] 
for the metallicity relations. For the Teff, we estimate a 
systematic uncertainty of 50 K. 

Using these methods, and combining the measurement 
and systematic uncertainties in quadrature, we mea- 
sure [Fe/H] -1-0.28 ± 0.14, [M/H] = -^0.20 ± 0.10 and 
Tefi = 3815 ± 88 K for KOI-254 using the iC-band analy- 
sis. Similar parameters appear in Paper 1, which uses the 
same data. The parameters are consistent with the re- 
sults from the photometric analysis within the estimated 
uncertainties. 

4. LIGHT CURVE ANALYSIS 
4.1. Joint Analysis and Planet Properties 

We simultaneously analyzed a total of 49 Kepler tran- 
sit light curves, a single Z-band light curve from the 
Nickel telescope at Lick Observatory, and our 14 RV 
observations acquired with Keck/HIRES. The Nickel 
and HIRES observation timestamps were converted to 
BJDc/TC to match Kepler MAST data using the tech- 
niques of Eastman et al. (2010). 

We fitted the Kepler and Nickel light curves using ver- 
sion 3.01 of the Transit Analysis Package (Gazak et al. 
2011), which uses the analytic eclipse model of Mandel 
& Agol (2002). For the Kepler transits we resampled the 
model to a cadence of 60 seconds before rebinning to the 
29.4 minute observing cadence to account for long inte- 
gration light curve distortions (Kipping 2010). We deter- 
mined the best-fitting parameters and their uncertainties 
using the same Metropolis-Hastings implementation de- 
scribed in § 3, with which we employ a Daubechies fourth 
order wavelet decomposition likelihood function (Carter 
& Winn 2009). Wavelet decomposition techniques pro- 
vide increased confidence in derived MCMC uncertainties 
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Fig. 7. — The Kepler (upper, blue) and Nickel (red, lower) light 
curves, phased at the photometric period. The Nickel light curve 
has been offset artificially for clarity. The best-fitting light curve 
models are shown for each data set (sec § 4), and the residuals are 
shown beneath each light curve. 



over the traditional likelihood by allowing parameters 
which measure photometric scatter (uncorrelated gaus- 
sian aw, and l/f correlated red Cr) to evolve as free 
parameters. The technique recovers the likelihood in 
the case where ar = and a^, is locked at a value charac- 
teristic to the observed data. For the RV data we fitted 
a Keplerian model using the partially linearized scheme 
of Wright & Howard (2009). 

Of the fifteen parameters in this technique, thirteen 
vary freely within our MCMC analysis: the period P, In- 
clination i, the scaled semimajor axis an, the radius ratio 
Rp/Rs, times of mid-transit Ttr, eccentricity e, argument 
of periastron uj, a^, fr, radial velocity amplitude K, the 
systemic velocity offset 7, and two parameters to account 
for global linear trends in the data normalization. The 
remaining two limb- darkening coefficients evolve under 
normal priors. For the Nickel Z-band data we adopted 
from Claret (2004): /if = 0.353±0.35, /i2 = 0.255±0.025. 
For the Kepler data we used the coefficients listed by Sing 
(2010): Ml = 0.521 ± 0.056, and ^12 = 0.225 ± 0.052. It is 
important to note that our joint fitting procedure allowed 
uncertainties in the orbital eccentricity to propagate into 
the determination of the Keplerian orbit parameters and 
the scaled semimajor axis an- 

We ran 40 independent MCMC chains each with 5 x 
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10^ links for a total of 1.4 x 10^ total inference links 
after removing the burn-in portion of the chains. Wc 
test for and find good convergence using the Gelman- 
Rubin statistic (Gelman & Rubin 1992). For the Kepler 
light curves we find typical values of = 0.00049 and 
CTr = 0.0018. For the Nickel light curve we find = 
0.0064 and cr^ = 0.024. 

Based on our light curve analysis, together with the 
stellar properties in § 3.7, we measure a planetary mass 
of Mp = 0.505 ± 0.090 Mjup, and a radius Rp = 
0.96 ± 0.19 Rjup. The complete list of stellar and planet 
parameters and their uncertainties is given in Table 4. 

We also tested for achromatic transit depths, 5, by 
fitting the Kepler and Nickel light curves separately. We 
found (^Kcpicr = 0.179 ±0.002 and Smckd = 0.183 ±0.016. 

4.2. Searching for Transit Timing Variations 

To measure the individual transit mid-times we fixed 
all of the global parameters [Rp/Ri,, a/Ri,, P, i and the 
limb- darkening coefficients) and fitted each transit event 
separately using the MCMC algorithm described in § 4. 
Table 3 lists the time at the mid-point of each transit 
Tmid; the difference between the measured values and 
those predicted by a linear ephemeris; and the formal 
measurement uncertainties, which are typically of order 
2 minutes. We see no statistically significant timing vari- 
ations. 

4.3. Limits on the System Age 

Measuring the ages of M-type stars is notoriously diffi- 
cult, except for the rare cases when stars are in clusters or 
associations, or show indications of extreme youth. One 
way of estimating the age of late-type field stars is gy- 
rochronology (Barnes 2010). Using this idea, age-period 
calibrations have been produced for G and K type dwarfs, 
(Mamajek & Hillenbrand 2008), but relationships for M 
dwarfs are still uncertain. 

The 122 day Kepler light curve shows clear evidence 
of rotational modulation, with 7 complete cycles shown 
in Figure 1. From these data we estimate a rotational 
period of 15.8 ± 0.2 days for KOI-254. This period can 
be directly compared with the rotational period sequence 
of the Hyades cluster (Delorme et al. 2011). From their 
Figure 15, and using our estimate oiV — K = 4.0±0.1 for 
KOI 254, it can be seen that stars in the 625 Myr Hyades 
cluster (Ferryman et al. 1998) with equivalent colors are 
rotating with a comparable average period of 14 days. 
This suggests that the system is relatively young, similar 
in age to the Hyades. 

We have generated a more detailed estimate of KOI- 
254's age using the gyrochronology relations derived by 
Barnes (2010). These relations express a star's cur- 
rent day rotation period as a fiinction of its convective 
turnover timescale and rotation period on the ZAMS. 
To calibrate our estimate, we adopted the relationship 
between a star's mass and convective turnover time as 
tabulated by Barnes & Kim (2010), and calculated three 
age estimates assuming the ZAMS rotation periods re- 
quired to reproduce the spread of rotation rates observed 
for ?a0.6 Mq stars in the 600 Myr Praesepe open cluster 
(Agiieros et al. 2011). Assuming a ZAMS rotation pe- 
riod of 2.81 days, as required to reproduce the median 
rotation period for O.6M0 Praesepe members, the Barnes 
(2010) gyrochrone relations predict an age of f«780 Myrs 



from the 15.8 day rotation period of KOI-254. However, 
a considerable range of periods are observed for Praesepe 
members in the mass range: adopting the ZAMS rota- 
tion periods required to reproduce the 10th and 90th 
percentile rotation periods observed for 0.6 Mq stars in 
Praesepe indicates that KOI-254's age could plausibly be 
as low as 380 Myrs or as old as 1.5 Gyrs. 

This estimate relies on the assumption that the ro- 
tational evolution of KOI 254 hasn't been grossly af- 
fected by tidal effects due to its hot Jupiter (e.g. Lanza 
2010). Prom its measured proper motion of 20.9 ± 
3.2 mas/yr (Monet et al. 2003) and our derived distance 
of 333 ± 33 pc, the system's tangential velocity is a mod- 
erate 34 ± 6 km/s. While deriving kinematic ages for 
individual star systems is unreliable, the space motion of 
KOI-254 indicates that it's age is consistent with 0.5 Gyr. 

5. SUMMARY AND DISCUSSION 

We report the detection of a short-period, Jupiter- 
mass planet that transits an early M-type dwarf star. 
The host star has mass M^, = 0.59 M0 and a radius 
Ri, = 0.55 R0, which we measure using a new method 
that draws upon broadband optical and NIR photometry, 
together with various empirical calibrations of M dwarf 
properties. Using our estimated stellar parameters, we 
find that the planet has a period of 2.455239 days, a 
mass Mp = 0.505 Atjup, and a radius Rp ~ 0.96 Rjup- 
At the planet's semimajor axis of 0.030 AU and given the 
star's estimated effective temperature of Toff = 3680 K, 
the planet has an equilibrium temperature of 1000 K (see 
Table 4 for the full list of parameters and uncertainties). 

Based on the rotation of the KOI-254, we estimate an 
age of roughly 0.5 Gyr. Examination of the tabulated 
theoretical planetary radii computed by Fortney et al. 
(2007) for an age of 0.3 Gyr and stellar insolation equal 
to that of KOT254 b shows that our measured planet ra- 
dius is consistent with model expectations, but only for 
a core mass of 50 M0. For the 1 Gyr models a 25 M^ 
core is necessary to reproduce the observed radius. This 
is a remarkably large core mass (heavy element content) 
and conforms well with the observed correlation between 
planetary heavy element content and host star metal- 
licity (Burrows et al. 2007; Torres et al. 2008; Miller & 
Fortney 2011). This in turn agrees with observations to 
date that suggest that the formation of a Jupiter-mass 
planet around a low-mass star requires high metallicity 
(Johnson & Apps 2009; Schlaufman & Laughlin 2010; 
Rojas-Ayala et al. 2010). 

As a Jupiter-mass planet with a period less than 
10 days, KOI-254 b is the sole example of a hot Jupiter 
orbiting an M dwarf. More than 300 stars have been 
monitored in various radial velocity surveys and a planet 
with P < 30 days has not been found around an M-type 
dwarf, despite the large expected Doppler signal. Indeed, 
a planet such as KOI-254 b would be readily detectable 
in any of the past and present M dwarf Doppler surveys 
conducted over the past 15 years (see e.g. Endl et al. 
2003; Johnson et al. 2010b; Bonfils et al. 2011). Given 
the difficulties associated with photometrically monitor- 
ing a large sample of M dwarfs, and the paucity of giant 
planets around low-mass stars, this fortuitous Kepler de- 
tection may serve as our lone example of a hot Jupiter 
around an M dwarf for some time to come. 
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Fig. 8. — Color-color diagram for low-mass stars in the Solar neighborhood. The "main sequence" is approximated by a constant for 
V — K < 5.5, and a polynomial for redder values (solid line). Metal-rich M dwarfs known to have [Fe/H] > -fO.25 as based on their their 
FGK dwarf companions are highlighted above the main sequence (orange circles). Metal-poor stars with [Fe/H] < —0.25 are shown below 
the main sequence (green circles), indicating a relationship between metallicity and the distance a star lies away from the main sequence. 



APPENDIX 

A. A NEW COLOR-BASED METALLICITY CALIBRATION 

Inspection of the locations of low-mass dwarfs in the (J — K)-{V — K) color-color diagram reveals a dramatic increase 
in the width of the main sequence beyond V — K k, A, corresponding to the onset of strong molecular band heads 
characteristic of the M spectral type (Figure 8). Several M dwarfs known to be metal- rich based on the [Fe/H] of 
their more massive binary companions lie along the upper envelope of the scatter in J — K . Similarly, several known 
metal-poor dwarfs reside in the lower portion of the diagram, suggesting that the increased scatter va. J — K color is 
due to the effects of metallicity. 

This color-metallicity effect has been noted previously by Leggett (1992) and Lepine & Shara (2005). The reason 
for the effect is most likely due to changes in continuous opacities due to molecular species such as Il2^, H^, He", 
H+ and C~ (AUard & Hauschildt 1995). Additionally, sources of line opacity can be seen in the spectral standards 
compiled by Rayner et al. (2009). The J-band (« 1.2-1.35 /im) spectra of M dwarfs exhibit deep potassium (K) and 
iron hydride (FeH) absorption features, along with dozens of shallower metal lines such as Na, Mg, Fe and Si. On the 
other hand, the iiTs-band (« 2.1-2.35 /xm) is relatively featureless, with a few relatively shallow Na and Ca lines as 
the only prominent absorption features. Thus, higher stellar metallicity preferentially suppresses J-band flux, causing 
the star's J — K color to become redder. 

To calibrate the relationship between J — K color and metallicity we first fitted a two-part function to the color-color 
diagram shown in Figure 8. This main sequence is approximated well by a constant {J — K)q — 0.835 for 3.8 < V — K < 
5.5. FotV -K > 5.5 we fitted a polynomial (J - K)o = X), a^iy - where {a} = {1.637, -0.2910, 0.02557}. We 
assume the main sequence in this color-color plane is an isometallicity contour with a value equal to the mean [Fe/H] 
of the Solar neighborhood. Based on an 18-pc volume-limited sample of stars in the Spectroscopic Properties of Cool 
Stars catalog (Valenti & Fischer 2005), Johnson & Apps (2009) measure a mean metallicity of [Fe/H] —0.05. 

Next, we follow the methodology of Bonfils et al. (2005) and use a calibration sample of M dwarfs with widely 
separated, FGK-type, common-proper-motion companions. The metallicity of the M-type companion in each pair is 
assumed to be the same as that of the more massive companion, for which spectroscopic parameters can be estimated 
accurately. We gathered our collection of 30 calibration stars from the SPOCS catalog and elsewhere in the literature. 

Figure 9 shows the metallicity of our calibration stars plotted against A{J — K) = (J — i<')„ioas~ {J ~ K)q. We fitted 
a linear relationship holding the offset fixed such that [Fe/H] = —0.05 at A(J — K) = 0, and found an acceptable fit 
of the form 



[Fe/H] = -0.050-f 3.520A(J-K-) 

= 2.872 + 3.520(7-/^) [ ±0.15 dex] (Al) 

valid for -0.1 < A( J - K) < 0.1 (A2) 

The residuals to our fit have a rms = 0.15 dex, which we adopt as the measurement uncertainty. It is important to 
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Fig. 9. — the metallicity of our calibration stars plotted against the distance of a star from the isometallicity contour shown in Figure 8, 
parametrized by A( J — K) = (J — -ft')moas — (J — K)q. The linear relationship (red solid line) is given by Eqn. A2. The dashed lines denote 
the A(J — K) limits between which our calibration is valid, roughly corresponding to —0.5 < [Fe/H] < +0.5. 

note that our calibration sample only spans —0.1 < A(J — K) < +0.1, corresponding to —0.5 < [Fe/H] < +0.5. 

We incorporate this calibration along with several other constraints in the analysis presented in § 3 to provide a 
jointly constrained calibration for stellar mass, metallicity and distance. To be consistent with our joint analysis, it is 
necessary to invert Eqn. A2 to give 

( J - if ) = 0.816 + 0.284[Fe/H] (A3) 

It is important to note that this simplified version of Eqn. A2 is only valid for our specific analysis of KOI-254, which 
hasV-K < 5.5. 
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TABLE 3 

KOI-254 Transit Mid-times and Ephemeris Residuals 





(BJD-2450000.0) 


Tmid - Ephemeris 


Telescope 


54964. 


5368 ± 0.0015 


-0.00048 ± 0.0018 


K 


54966. 


rinooo 1 n r\nr\oi 

,99z28 ± U.UUUoi 


-U.UUU24 ± U.UUio 


K 


54969. 


,44Dyo ± U.UUUo4 


-U.UUU/D ± U.UUlo 


Iv 


54971. 


.yuoUo ± u.uuuyi 


U.UUUUo9 ± U.UUi4 


T/" 

Iv 


54974. 


.ooooo ± U.UUUyo 


U.UUUlO ± U.UU14 


T/" 

iv 


54976. 


.olo57 ± (J.UUUyi 


U.UUU14 ± U.UU14 


K 


54979. 


.zbooo ± U.UUUoo 


n nnnQQ _L n nni q 
-U.UUUoo ± U.UUlo 


TV" 

rv 


54981. 


.72302 ± U.0U096 


n nnnoT i n nm A 

-0.00087 ± 0.0014 


K 


54984. 


.17923 ± U. 00086 


0.00011 ± 0.0013 


K 


54986. 


.6343 ± 0.0015 


-0.000025 ± 0.0018 


K 


54989. 


.08909 ± 0.00091 


-0.00049 ± 0.0014 


K 


54991. 


,54493 ± 0.00068 


0.00012 ± 0.0012 


K 


54993. 


.9996 ± U.UUlU 


-U.UUU40 ± U.UU14 


K 


54996. 


,4o40D ± U.UUU9D 


n nnn'7n _L n nm /i 
-U.UUU/U ± U.UU14 


Iv 


55003. 


.8211 ± 0.0010 


n nnm c i n nm /i 

0.00015 ± 0.0014 


K 


55006. 


.27645 ± 0.00085 


n nnnT7 i n nm o 

0.00027 ± 0.0013 


K 


55008. 


.73088 ± 0.00085 


n nnncT /i i n nm o 

-0.00054 ± 0.0013 


K 


55011. 


.18647 ± 0.00081 


n nnm t i n nm o 

-0.00017 ± 0.0013 


K 


55013. 


,d4zU ± U.UUiU 


n nnm o i n nm /i 

U.UUUi2 ± U.UUi4 


K 


55018. 


.OOZiU ± U.UUU9o 


n nnnoQ _L n nm a 
-U.UUU20 ± U.UU14 


Ty 
iv 


55021. 


.UUo4 ± U.UUiU 


n nnnoo _L n nm a 
U.UUUoo ± U.UU14 


yy 

IV 


55023. 


.4dx4z lb U.UUUo4 


-U.UUUoo ± U.UUlo 


TV- 

rv 


55025. 


.91715 ± U.UUUTo 


n nnnc? _!_ n nm o 
-U.UUUo7 ± U.UU12 


TV- 

rv 


55028. 


.37311 ± 0.00091 


n nnm d i n nm o 

-0.00014 ± 0.0013 


K 


55030. 


.82892 ± 0.00083 


0.00044 ± 0.0013 


K 


55033. 


,2832 ± U.0U2U 


-U.UUU51 ± U.UU22 


K 


55035. 


. / o9o ± U.UUil 


n nnni^n _L n nm c; 
U.UUUu9 ± U.UUlo 


iv 


55038. 


.19468 ± 0.00084 


n nnncro i n nm o 

0.00052 ± 0.0013 


K 


55040. 


.04897 ± U.UUUo2 


n nnn /I o i n nm o 

-U.UUU42 ± U.UU12 


K 


55043. 


.iU4iO ± U.UUUoo 


n nnm q _l n nm q 
U.UUUio ± U.UUlo 


Iv 


55045. 


.0090D ± U.UUU9i 


n nnnnn _L n nm a 
-U.UUU29 ± U.UU14 


is. 


55048. 


.UioUo ± U.UUUoi 


n nnnnc; _L n nm q 
U.UUU90 ± U.UUlo 


Iv 


55050. 


.47094 ± 0.00093 


n nnni^o i n nm /i 

0.00063 ± 0.0014 


K 


55052. 


.92573 ± 0.00088 


n nnnnn i n nm o 

0.00020 ± 0.0013 


K 


55055. 


.3827 ± 0.0011 


n nm n i n nm c 

0.0019 ± 0.0015 


K 


55057. 


.83570 ± 0.00070 


n nnnon i n nm o 

-0.00030 ± 0.0012 


K 


55060. 


.29115 ± 0.00093 


n nnnnTo i n nm a 

-0.000073 ± 0.0014 


K 


55062. 


.74626 ± 0.00089 


n nnm n i n nm o 

-U.UUU19 ± U.UUlo 


K 


55065. 


,20262 ± 0.0008; 


n r\r\r\CiA _L n nm q 
U.UUU94 ± U.UUlo 


IV 


55067. 


,65679 ± 0.00091 


-0.00012 ± 0.0014 


K 


55070. 


,11231 ± 0.00071 


0.00017 lb 0.0012 


K 


55072. 


,56725 d= 0.00085 


-0.00012 ± 0.0013 


K 


55075. 


,02311 ± 0.00089 


0.00051 lb 0.0013 


K 


55077. 


,47754 ± 0.00097 


-0.00029 lb 0.0014 


K 


55079. 


,9343 ± 0.0010 


0.0012 lb 0.0014 


K 


55082. 


,38913 ± 0.00094 


0.00084 lb 0.0014 


K 


55084. 


,84365 ± 0.00086 


0.00013 zb 0.0013 


K 


55087. 


,29916 ± 0.00084 


0.00041 lb 0.0013 


K 


55089. 


,75411 ± 0.00099 


0.00014 lb 0.0014 


K 


55742. 


,84449 ± 0.0027 


-0.00047 lb 0.0029 


N 



Note. — K — Kepler, N — Nickel Z-band 
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TABLE 4 
System Parameters for KOI-254 



Parameter 


Value 


68.3% Confidence Comment 






Interval 



Transit Parameters 
Orbital Period, P [days] 
Radius Ratio, (Rp/Ri,) 
Transit Depth, [Rp / R^,)'^ 
Scaled semimajor axis, a/R^ 
Orbit inclination, i [deg] 
Transit impact parameter, b 
Rotation Period, Prot [days] 

Other Orbital Parameters 
Eccentricity 

Argument of Periastron u) [degrees] 
Velocity semiamplitude i^* [m s~^] 

Stellar Parameters 

M„ [Mq] 

R* [R@] 
p* [pe] 
logff* [cgs] 
[Fc/H] 
[Fe/H]NiR 
[M/H]nir 
Distance [pc] 

Teff.NIR [K] 

Planet Parameters 
Mp [Afj„„] 
Rp [Rjup] 

Mean planet density, pp [g cm~^] 
logffp [cgs] 

Equlibrium Temperature Tefj (il*/a)i/2 [K] 



2.455239 ±0.000004 A 

0.179 ±0.002 A 

0.03204 ±0.00072 A 

10.6 ±1 A 

87.0 ±0.7 A 

0.6 ±0.2 A 

15.8 ±0.2 A 

0.11 t°o lg A, C 

230 ±68 C 

110 ±10 C 

0.59 ±0.06 D 

0.55 ±0.11 D 

4 ±2 A 

4.7 ±0.2 B 

+0.13 ±0.13 D 

+0.28 ±0.14 F 

+0.20 ±0.10 F 

333 ±33 D 

3820 ±90 G 

0.505 ±0.090 B,C 

0.96 ±0.11 B 

0.8 ±0.5 B,C 

3.1 ±0.1 A 

1000 ±20 D 



Note. — Note. — (A) Determined from the light curves. (B) Based on group A pa- 
rameters supplemented by the photometric stellar mass determination described in § 3. 
(C) Based on our analysis of the Keck/HIRES RV measurements. (D) Based on our pho- 
tometric mass and radius determinations described in § 3.3. (E) Based on photometric 
metallicity calibrations descrbed in § 3 and the Appendix. (F) Based on the analysis of 
our TripleSpec NIR spectra using the methodology of Rojas-Ayala (2011). 



